Weak localization in ferromagnets with spin-orbit interaction 
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Weak localization corrections to conductivity of ferromagnetic systems are studied theoretically 
in the case when spin-orbit interaction plays a significant role. Two cases are analyzed in detail: 
(i) the case when the spin-orbit interaction is due to scattering from impurities, and (ii) the case 
when the spin-orbit interaction results from reduced dimensionality of the system and is of the 
Bychkov-Rashba type. Results of the analysis show that the localization corrections to conductivity 
of ferromagnetic metals lead to a negative magnetoresistance - also in the presence of the spin-orbit 
scattering. Positive magnetoresistance due to weak antilocalization, typical of nonmagnetic systems, 
does not occur in ferromagnetic systems. In the case of two-dimensional ferromagnets, the quantum 
corrections depend on the magnetization orientation with respect to the plane of the system. 
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PACS numbers: 72.25.-b, 72.15.Rn, 73.20.Fz 
I. INTRODUCTION 



Owing to the giant magnetoresistance effect discov- 
ered in artificially layered structures, Em transport prop- 
erties of low-dimensional magnetic systems were exten- 
sively studied in the past decade. The huge interest was 
stimulated by applications of the effect in magnetoelec- 
tronics, particularly in read/write heads, field sensors, 
random access memory elements, and othersJ3 Since the 
effect exists also at room temperatures (which is impor- 
tant for applications), there was only a little interest in 
the low temperature regime, where quantum corrections 
to conductivity may play a certain role. This regime, 
however, may be important in view of possible applica- 
tions of metallic a*id/or semiconducting macsietic sys- 
tems in spintronicsu and quantum computings 

The quantum corrections to electrical conductivity, re- 
lated to scattering of electrons from impurities in non- 
magnetic metals and doped semiconductors, were exten- 
sively studied in the past two decades.Bij However, the 
problem of quantum corrections in ferromagnetic met- 
als is still unexplored. Only a few works on this sub- 
ject can be found in the relevant literature. These in- 
clude two theoretical worksEIM and a few reports on 
experimentsjlailll which prove the existence and impor- 
tance of the quantum corrections related to both weak 
localization and electron-electron interaction effects. The 
theoretical description, however, was restricted to the ef- 
fects of localization on the spin-density ilnctuations in the 
vicinity of the ferromagnetic transition!!!] and to electron- 
electreji interaction effects in spin dependent quantum 
wellsH 

It is well known that the quantum corrections due to 
weak localization in nonmagnetic systems are suppressed 
by a sufficiently large magnetic induction B. One may 



then expect a similar suppression of weak localization 
by an internal magnetic induction Bi nt in ferromagnets. 
But this point is still not definitely clear from the exper- 
imental point of view, at least for some kinds of ferro- 
magnetic materials. It is then reasonable to assume that 
the internal magnetic induction existing inside the ferro- 
magnets may reduce the localization corrections instead 
of destroying them totally. Very likely, one can expect 
only a slight effect of Bi nt in the case of newly .xbsveloped 
magnetic semiconductors like GaMnAs alloysEj 

Spin-orbit (SO) scattering from paramagnetic impu- 
rities in nonmagnetic metals is known to have a signifi- 
cant influence on the quantum corrections. It can reverse 
the sign of the localization correction (so-called weak an- 
tilocalization effect), which results in_a .positive magne- 
toresistance at weak magnetic fields£jcJ However, SO 
interaction may also result from other sources—like for 
example the DresselhausH or Bychkov-RashbanJ terms 
in the relevant Hamiltonian. These terms are related to 
the lack of inversion symmetry in certain crystals or to 
reduced dimensionality in quantum-confined structures, 
respectively. In the context of weak localization theory, 
this type of spin-orbit interaction has been^-studied by 
A. G. Aronov and Yu. B. Lyanda-GellerHta 

Recently, there is a large interest in SO interaction in 
magnetically ordered materials. First, the SO scattering 
is believed to be responsible for the anomalous Hall effect 
in ferromagnets. Ol3 Second, the SO interaction is one 
of the main interactions which determine the spin diffu- 
sion length. The latter quantity plays a crucial role in 
the dependence of the giant magnetoresistance effect on 
the sublayer thicknesses., when the current flows perpen- 
dicularly to the films-BEl 

In this paper we study the localization corrections to 
conductivity of ferromagnetic systems in the presence 



1 



of SO scattering from defects and also in the presence 
of the Bychkov-Rashba term. It is well known that in 
nonmagnetic materials the spin-orbit scattering is cru- 
cial for the localization correction. As it was already 
mentioned above, the SO scattering leads to weak an- 
tilocalization, im~ .to positive magnetoresistance at small 
magnetic fields .li^iEll The situation in ferromagnetic met- 
als, however, is significantly different. We show that the 
processes, leading to weak antilocalization in nonmag- 
netic systems, are totally suppressed in ferromagnets, so 
in the presence of SO interaction we have only a negative 
magnetoresistance. 

In Sections 2 and 3 we derive the formula for the 
Cooperon and quantum correction to conductivity in 
three-dimensional (3D) ferromagnets. Two dimensional 
(2D) ferromagnets are considered in Section 4, while 
quantum wells and possible crossover from 2D to 3D case 
are discussed in Section 5. The influence of SO interac- 
tion in the form of Bychkov-Rashba term is studied in 
Section 6. In Section 7 we discuss the effect of internal 
magnetization on the localization corrections. Finally, 
conclusions and summary are in Section 8. 



II. COOPERON IN A 3D FERROMAGNET 

We consider the following model Hamiltonian of a fer- 
romagnet with SO scattering: 



H = d 3 r^ (r) 



V 2 

-— -Ma z +V(r) 
2m 



1><?) > (1) 



where the axis z is assumed to be along the magnetiza- 
tion M, -0 is a spinor field with the components if}-\ and 
i{>l, and we put h = 1. In the presence of a magnetic in- 
duction B = rot A, the gradient operator V is replaced 
by V — i eA/c. Note that the exchange term —Ma z acts 
only on the spin and has therefore no direct effect on the 
orbital motion of the electrons. 

The random potential V(r) of impurities consists of 
two statistically independent components - the compo- 
nent independent of the electron spin, described by the 
random potential Vo(r), and the spin-orbit component 
V so (r). Matrix elements of the latter component have 
the form 



(V r so )k Q ,k'/3 = i Vi (k x k') • a a p 



(2) 



for the transitions (k, a) — > (k',/3), where V\ is a con- 
stant, k and k' are the initial and final electron wavevec- 
tors, a and (3 describe the corresponding spin states, and 
g = {&x, Py,&z) are the Pauli matrices. 

The key conmonent of the weak localization theory is 
the Cooperon,Qia which can be presented by a ladder 
in the particle-particle channel with two propagators de- 
scribing electrons with almost vanishing total momentum 
and with very close energy parameters (so-called Cooper 
channel). In the case of ferromagnets, and as long as 



M 3> t^T 1 , rj~ 1 , where v\ and T[ are the momentum re- 
laxation times of the spin up and down electrons at the 
Fermi surface, this channel does not include ladder ele- 
ments with Green functions corresponding to the oppo- 
site spin orientations. Indeed, y&ing the standard method 
of calculation of the Cooperon,EH3 one should evaluate the 
integral 
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d 3 k 

(2^)3 



G?(e,k)G£(e,-k), 



where G^' A (e, k) are the retarded and advanced Green 
functions of an electron with spin a =T,|- This gives 

n T{ /n TT ~ i/Mr T < i, and n n /n u ~ i/M n « 1, 

which corresponds to vanishingly small contribution of 
the singlet channel. This result can also be understood 
as a suppression of singlet pairs by the exchange field. 

Validity of this approach is confirmed by the following 
estimations. If we assume parameters typical for pure 



Fe,E3 i.e., M = 2.5 eV and r T ~ tj. ~ 5 x KT 13 s, we ob- 
tain (Mrj,|) _1 ~ 5 x 10~ 4 . In dirty Fc, this value would 
be increased by one or two orders of magnitude, but will 
still remain small as compare to unity. The exclusion of 
the Cooperon in the singlet channel is the crucial point 
of our description, which leads to the absence of weak an- 
tilocalization in ferromagnets. In the following, we will 
omit the spin index in the inequalities M t7 ,t7 , 



and will write simply M ^> r 



a, k 
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FIG. 1. The bare scattering amplitude. 

In the case of weak scattering potential and upon aver- 
aging over the random field V(r), one finds the following 
bare scattering amplitude (see Fig. 1): 

F Vy<5 = aS a55p~f - b (o-% S <7fa + vls^tf) - taPiS , (3) 

where i a/ 3 7 a is a diagonal matrix in the space of (|T)- (Tl 
),(IT),ai) states, 

t = diag(di, d 3 , d 3 , d 2 ) , (4) 

and we introduced the following definitions: 
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a = Ni Vj = 



2tti/ t t ot 271-^Toi 



b = Ni Vl X Q fc| T k 2 Fl , 



(5) 



(6) 
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di = Ni V 2 A 4 T = 



d 2 = Ni V 2 A k 4 Fl 



d 3 = -Ni V 2 A k 2 n k 2 Fl , 



with An defined as 



A = (n k x n k , 



(7) 



(8) 



(9) 



(10) 



In the above equations Ni is the concentration of scat- 
tering centers, and are the densities of states at the 
Fermi level for spin-up (majority) and spin-down (minor- 
ity) electrons, while kp\ and fcpj are the Fermi wavevec- 
tors for spin-up and spin-down electrons. The relaxation 
times Tof an( l r oj.i defined by Eq. (5), are the momentum 
relaxation times in the absence of SO scattering, whereas 
the relaxation times r* + and r* r, defined by Eqs. (7) 
and (8), are due to the SO scattering. The averaging in 
Eq. (10) is over all orientations of the unit vectors n k 
and n k /, where n k = k/fc and n k / = k'/fc'. 

Using Eq. (3), we find the following bare scattering 
amplitudes for the spin-up and spin-down electrons: 



TTTT 



= a — d\ = 



1 



1 



27TI/|To-f 27TZ/JT 



(11) 



Summing up the ladder diagrams contributing to the 
renormalized scattering amplitude r-|-(a;, q) = T^^(lu, q) 
in the Cooper channel with small transferred energy u 
and momentum q Qu (see Fig. 2), one obtains the equa- 
tion 



where 



r T (w, g ) = r l T ) + r l T ) n T (o;,( 7 )r T ( w , (Z ) , (13) 



n T (w,«) = J 7^3 Gf (w,k + q) Gf (0, -k) 



27r^-fT| (l — D^q 2 r^ + iujrA 



(14) 



Here, £)f = ^ w^Tf is the diffusion coefficient, vp-\ is the 
Fermi velocity, and tj is the momentum relaxation time 
of spin up electrons. Equation (14) was derived in the 
diffusion limit, i.e., when uj <C 1/r-f and g <C (Z?|T|) -1 / 2 . 
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l T sol 



(12) 



2^1 r oi 

One should note at this point that the bare elements 
and r°-|.^-|. do not contribute to the Cooper-channel 
diagrams in the case of ferromagnetic systems as long as 
M >> 1/t, as we already stated before. Apart from this, 
in a 3D case T° does not contain the components r° 



TTli 



and r^fT' which vanish due to the rotational symmetry 
in the x-y plane, as can also be concluded from Eq. (3). 



FIG. 3. Self energy diagram. 

The relaxation time tj can be found by calculating the 
self energy, presented by the diagrams of Fig. 3. The self 
energy contains non-zero spin-flip vertices of the singlet 
type, as shown in Fig. 4. After calculating the imaginary 
part of the self energy, we find 



— = 2ttv^ I a + d\ 

n 



2n 



(15) 
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FIG. 2. Ladder diagram for the Cooperon in 3D case. 



FIG. 4. Spin-flip vertex. 
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Using Eqs. (5) to (14), we finally obtain the following 
formula for the renormalized scattering amplitude in the 
Cooper channel (Cooperon): 

1 a + di + (2vjv^)b 



r T (w,q) 

3 effec 
is introduced, 



(16) 



-iu; + D^q 2 + l/f soT 
where the effective SO relaxation time of Cooperon, f so ^ , 



J_ = 2_ d 1 + {v l /v 1 )b 
f so f T| a — di 



(17) 



Let us now define rf ^ and as 1/t£ ^ — and 
1/t^ i = 27rz^|&. In the limit of weak SO scattering, 
r so| r T ' 0110 ma y then write 



r T (w,g) 

where 



27r^ T r-p -iw + D^q 2 + l/f so t + 1/r, 



= 2 



1 



1 



T so \ 



(18) 



(19) 



and a phase relaxation time t v -[, related to some in- 
elastic scattering processes of electrons pB is added into 
Eq. (18). 

The analogous formulae can also be derived for the 
spin-down Cooperon Ti(u,q) = T ml (uj,q). This for- 
mula can be obtained from Eqs. (18) and (19) by invert- 
ing direction of the arrows. 



III. CONDUCTIVITY AND 
MAGNETO-CONDUCTIVITY IN 3D CASE 

The localization correction Act ter-fhe static conductiv- 
ity is determined by the diagramsB'Q shown in Fig. 5. 





FIG. 5. Localization corrections to conductivity. 

Upon calculating the corresponding contributions, one 
finds the following formula for Act: 



Act = -7 (2^ T r T 2 ^T / ^T 7 (0, q )+2^ l rfD l J ^2^(0, 



a) 



(20) 



which is a straightforward generalization of the localization correction in a nonmagnetic case. Using Eqs. (18) and 
(20), we obtain 



Act = const - 



4tt 2 



D 1 ' 2 \ T so-\ T v t 



1/2 



L) i \ T so[ 



1 



1/2 



(21) 



where the constant part is related to the contribution from theJargest momenta of Cooperon, q ~ {Dt) 1 I 2 , and can 
not be calculated exactly within the diffusion approximation.Bi3 It can be estimated as 



const 



e 

~4^2 



(Lv T r 1/2 + (^nr 1/2 



Thus, the constant term is negative and, since tjj < f so , t v , it is larger in magnitude than the second term of Eq. (21). 
Therefore, the total correction (21) is negative. By decreasing f so and/or t v , we suppress the localization correction 
to conductivity. 

The magnetic induction (both external and internal) supresses the localization coraections. If the magnitude of the 
total magnetic induction is B, then, following the method developed by KawabatajEl we find 



Aa(B) 



eB 

27TC 



"■OT 

E 



dq 
2^ 



n n 2 



AeBDt 



1 



1 

1> T 



D iq 2 



4eBDi 



U+ 2 



1 

Tsol 



1 

T <pl 



(22) 
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where Ib = (c/eB) 1 / 2 is the magnetic length, c is the light velocity, and the sums over the Landau levels are cut off 

atnot(i) ( D Ul) T Ul))- r-i 

After eliminating the Aa(B = 0) part JEj we find a formula, which is a generalization of the Kawabata's low-field 
magnetoresistance to the ferromagnetic case, 



Aa(B) - Acr(0) 



16tt 2 Z, 



■E 

n=0 



(n + l/2 + * T )V 



— -2{n + l + 5^) 1 l 2 + 2{n + 5 ] ) 1 / 2 



(n+ l/2 + ^)V2 



-2{n+\ + 5 i ) 1 ' 2 + 2{n + 5 i ) 1 /' 



(23) 



where 



>TU) 



I 2 
'b 



4D 



TU) 



«°TU) 



(24) 



In accordance with Eq. (22), the magnetic induction su- 
presses the negative correction to the conductivity. Thus, 
the resulting sign of Eq. (23) is positive, and its magni- 
tude, Aa(B) — Acr(0), increases with increasing magnetic 
field. This means that one finds a negative magnetore- 
sistance, despite of the presence of the spin-orbit inter- 
action. The reason of this is the fact that we have ex- 
cluded the singlet Cooperon, which contributes to the 
localization correction with the opposite sign, and usu- 
ally gives rise to a positive magnetoresistance (weak an- 
tilocalizatkgLl in weak magnetic fields in nonmagnetic 
materials Enu 

The obtained result is in agreement with the results on 
magnetoconductivity of nonmagnetic metals in a mag- 
netic field, if both Zeeman splitting and spin-orbit scat- 
tering are taken into account.EllB Indeed, the exchange 
field of a ferromagnet enters the Hamiltonian, Eq. (1), 
like the Zeeman term. Thus, the strong exchange-field 
limit of a ferromagnet corresponds to the case of a large 
magnetic field in the Zeeman term. At these conditions, 
the effect of a magnetic field is associated with a negative 
magnetoresistance due to the suppression of the singlet 
Cooperon by magnetic field through the Zeeman split- 
ting. 



IV. TWO-DIMENSIONAL FERROMAGNETS 



In this section we consider a two-dimensional ferromag- 
net. In such a case there is no electron motion in the di- 
rection perpendicular to the plane, and consequently the 
electron wavevectors are in the plane of the ferromagnet. 
We consider first the case of in-plane magnetization, as 
shown schematically in Fig. 6. 



M 



FIG. 6. 2D ferromagnet with in-plane magnetization. 



The bare scattering amplitude has then the form 



where 



and 



2'KVTn 



27WT? 



(25) 
(26) 

(27) 



In a strictly 2D case and when both spin sub-bands are 
populated with electrons, the density of states is indepen- 
dent of the spin orientation, = = v. This follows 
from the model assumed in Eq. (2). In that case also 
the relaxation times in the absence and/or presence of 
the SO scattering are independent of the electron spin, 
Tot = T oi = T o and 7f ot = 7f ol = 7f . 

Contrary to the 3D case considered above, we have 
now to include the spin-flip processes in the Cooperon 
ladder. This is due to the fact that now we do not have 
the rotational symmetry in the x-y plane. According to 
Eqs. (25) to (27), the bare amplitudes are now, 



r? 



1 



27W7o 



1 TU1 - 1 UTT - 1 sf - 



1 



(28) 



(29) 
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The equations for the renormalized vertices can be writ- 
ten as two coupled ladders for (u>, q) =Tnn(uj,q) = 
T(uj,q) and r TTU (u;,g) = T lm (uj,q) = V sf (uj,q), as 
shown in Fig. 7. 





FIG. 7. Equation for the Cooperon in 2D ferromagnet with 
in-plane magnetization. 

From these equations we find 



4ttvt 2 -iu + Dq 2 + 1/t v ' 



(30) 



where D = \ (Dj + DjJ, and the diffusion constants 
and Di are defined as -DjQ) = 5 u Fjm T - The electron 
relaxation time r in Eq. (30) is independent of the elec- 
tron spin, Tf = T[ = t, and can be calculated in the same 
way as in the 3D case, which gives 



I 



I 

n 



i 

T 



1 

TO 



(31) 



It is worth to note that the spin-orbit scattering enters 
the Cooperon only through the one-particle relaxation 
time r, and has no influence on the pole of the Cooperon. 



FIG. 8. 2D ferromagnet with magnetization perpendicular 
to its plane. 

Consider now the case when the magnetization field 
M is perpendicular to the plane of the ferromagnet, as 
shown schematically in Fig. 8. The calculations are sim- 
ilar to those in the case of in-plane magnetization, so 
we write down only the results. The only bare scatter- 
ing amplitudes arc now and rj, which generally are 
different and have the form 



MU) 



l 



1 



(32) 



27WT0 — < soT(i) 

The Cooperons T^(u>,q) and T±(l}, q) can be written as 
1 1 



TU) 



27WT T(i) 



-iu + D Ui) q 2 + l/f soTa) + l/r vT(i) ' 

(33) 

where the relaxation times tj and are given by 



r tU) 



TO + 



sot a) 



and 



T s°TU) T so T (4) _ 2t T(4) 'soUD 

provided that t* ^^ 3> TjQ). By comparing the results 
(30) and (33), both obtained for a 2D case but for dif- 
ferent magnetic configurations, we see that the effect of 
SO interaction significantly depends on the orientation 
of the magnetization field M with respect to the plane of 
the ferromagnet. 

The quantum correction to the conductivity of a 2D 
ferromagnet takes the form 



(34) 



(35) 



I corresponding for- 
El Here, r^m and 



which is a direct generalization of tl 
mula in the 2D nonmagnetic case. 
t so |Q) are dehned, respectively, by (31) and t so |Q) = 
for the in-plane magnetization, and by Eqs. (34) and (35) 
for the case of perpendicular magnetization. 

The 2D localization correction, described by Eq. (36a), 
is negative since r < f so ,T v , and, in addition, we take 
here r <C t so ,t v . The latter inequality means that the 
momentum relaxation time of electrons, r, is mainly due 
to the potential scattering. 

We can also present an expression for the conductiv- 
ity in the case of nonzero magnetic induction B, per- 
pendicular to the plane, by generalizing the result for a 
nonmagnetic two-dimensional systemEJ 



i 



i 



2 T|0| 



1 



1 



?~soT a T 



r ¥>T a T 



The situation with a large number of occupied subbands 
is typical for metals, whereas the situation with only a 
few populated subbands is characteristic of semiconduc- 
tor quantum wells. The effect of SO scattering on the 
magnetoresistance of nonmagnetic materials at such con- 
ditions has been studied by Hikami et al.c3 

We consider here a ferromagnetic quantum well in the 
geometry shown in Fig. 6, i.e., when the magnetization 
field M is in the film plane. The scattering amplitude 
has then the following form 

r °af3jS = a Sasdp-y - b a% s ct| 7 - c a v aS a v ^ - t afjlS , (37) 

with a, b and ta^s defined by Eqs. (4) to (6), and 



c = Ni V{ (k x k') 



(38) 



n a i 



i 



Tsoiai r v iai 



(36b) 



= AeBD^i/c, and tp(x) is the Digamma 
which has the property ip(x) ~ ln(x) for 



where a* 
function! 
x > 1. 

The magnetic induction suppresses the negative cor- 
rection to the conductivity which leads to the negative 
magnetoresistance. It should be noted that in strongly 
2D case, the in-plane magnetic induction does not affect 
the localization correction to conductivity. The reason 
is that in the two-dimensional case, the flux of magnetic 
induction does not penetrate through any closed electron 
paths. Correspondingly, the in-plane induction does not 
break the interference of closed trajectories of electrons 
moving in opposite directions, responsible for the weak 
localization effect Era 



di = N t V{ (k x k'). 



r(T) 



Ni Vi (k x k') 



-(1) 



2 ™T T so T 



d 3 = -Ni vt (k x ko: 



t(T,1) 



(39) 



(40) 



(41) 



The averages in Eqs. (38) to (41) include averaging over 
discrete subbands due to quantization of the electron mo- 
tion along the axis x (normal to the film plane). In this 
notation, the expressions for bare vertices T9 and co- 
incide with Eqs. (11) and (12), whereas for the spin-flip 
vertex T° s j: we have 



V. QUANTUM WELLS 

In a quasi 2D case the electrons are confined within 
a quantum well and the number Np of 2D subbands 
populated with electrons is larger than one, Np > 1. 



1 sf 



27r ^T T foT 27W T T soT 



(42) 



The ladder equations for T^(ui, q), T^(uj,q) and 
T s f(to,q) have the form shown in Fig. 7. From these 
equations one finds the following solution for T-^(u), q): 



r T (w,q) = 



(ro / ) 2 n i ( w , q ) + ro i-r°iL>, 9 ) 



l - r° n T (w, g )j [i - r° n x ( u , 9 )j - (r° f y n T ( w , <?) n x ( u , q) 

In the case of weak SO interaction, the final expression for T^(u),q) takes the form 



r T (w,q) 



-4 1 



where 



Ar = 



2/^ oi + l/7f oT + l/r» T 



(43) 



(44) 



(45) 
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D 



D t = 



(2/7^ + l/rf oT + l/rf oT ) + ZJ X (2/r a ' oT + + 

2/t/ot + 2 A? o1 + (V^ot + V^t) (! + nM) 



(46) 



7"soT 



/ (^oT <) + (V^oT + VjW) (V^T + »U_ ("T ^ oi )) + 2/ (g T r^) 
2/^ oT + 2/t* A + (l/r* T + l/rJ oT ) (1 + ^/ n ) 



(47) 



and the spin-up relaxation time tj is 



1 



- = — + z^ + z^ + n7-- (48) 



tot 



■sol 



■sol 



The corresponding expression for Fj(cj,q) can be ob- 
tained from Eqs. (44) to (48) by changing the arrow direc- 
tion in these formulas. It is worth to note, that according 
to Eq. (46), the spin-orbit interaction renormalizes the 
diffusion coefficients. 

The correction to conductivity is described either by 
Eq. (21) or by Eq. (36), in dependence on the ratio of 

(~ \l/2 

L e ff — [D/(f~ 1 +T" 1 )) to the width of the quan- 
tum well L. For L e // < L we have effectively a 3D case 
described by Eq. (21), whereas for L e // 3> L one finds 
effectively 2D behavior, Eq. (36a). It is worth to note 
that due to a strong asymmetry between spin-up and 
spin-down states, a "mixed" case is possible, with a 3D 
correction for one spin orientation and a 2D correction 
for the opposite spin orientation. 



VI. SPIN-ORBIT INTERACTION OF THE 
BYCHKOV-RASHBA TYPE 



The case when the spin-orbit coupling enters the 
Hamiltonian also in the absence of scattering defects 
needs a special treatment. In the following we will con- 
sider a system described by the Hamiltonian including 
the Bychkov-Rashba interaction termnJ 



H = 



d 3 k 

(2tt) 3 



4 



2m 



Ma z + g (k y a z - k z a y ) 



(49) 

where the magnetization M is assumed to be in the plane 
of the system (the geometry of Fig. 6) and g = g(k) is 
the spin-orbit coupling parameter. This Hamiltonian can 
describe electrons in low-dimensional structures close to 
interfaces, or electrons in a quantum well with variable 
doping. 

The energy spectrum of the Hamiltonian (49) has two 
branches, 



e i (a) (k) = ^±[(M. 



gk y ) +g 2 k 2 z 



1/2 



(50) 



which are no longer pure spin-up and spin-down states, 
but correspond to spin-mixed states. The eigenfunctions 
corresponding to the eigenvalues (50) can be written as 



I lk) = (4M 2 +g 2 k 



H<?Mk T> + 2M|k 1)) 



(51) 



1 2k) = UM 2 



g 2 k 2 y 1/2 (2M\kV +l gk z \kl)) . (52) 



Due to the terms linear in k y and for M ^ 0, the en- 
ergy spectrum (50) is not symmetrical with respect to 
the k — + — k transformation, as can also be seen in Fig. 9, 
(a) and (b). 




2 4 
k z , 10 cm" ^ 




2 4 
k y , 10 cm" 



FIG. 9. The energy spectrum E(k) of a ferromagnet with 
the Bychkov-Rashba SO coupling, shown as a function of k z 
at k x = k v = (a) and as a function of k y at k x — k z — (b). 

In the absence of external electric field E, there are 
nonvanishing spin currents, associated with each branch 
of the spectrum Ei j a(k). However, it can be verified that 
the total charge current is zero if E = 0. 

In the basis of the spin-up and spin-down states, the 
electron Green function has the matrix form, 



Go (e,k) = 



£ — k 2 /2m + fx — Ma z + gk y a z — gk z a y 

{s - £ik + A* + ifi signs) (e - 62k +fi + iS signe) 

(53) 

where /1 is the chemical potential. Accordingly, the self 
energy is also a matrix in this basis. It is, however, more 
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convenient to consider the electron self energy S(e,k) 
in the basis of the eigenfunctions (51) and (52). The 
self energy is then diagonal and the imaginary parts of 
£i(0, kFi) and £2(0, &F2) give the relaxation times t\ 
and T2 in the presence of defects. 

Using Eqs. (51) and (52), one can calculate the matrix 
elements of the impurity potential Vb (r) , 



Ik' 



V 2k> 



2k' 



V {AM 2 +gg'k z k' z ] 



{AM 2 + g 2 k 2 z ) 1/2 {AM 2 + g' 2 k' 2 ) 



1/2 ' 
(54) 



V, 



lk, 2k' 



2k,lk' 



2iV M {gk z + g'k' z ) 



{AM 2 + g 2 k 2 ) 1/2 {AM 2 + g' 2 k' 2 ) 1/2 

(55) 

where g' — g{k'). 

To find the self energy one needs to calculate the dia- 
grams shown in Fig. 10. 



1 k 



1 k 



1 k 



1 k 



2 k 



1 k 



FIG. 10. 
interaction. 



Self energy in the case of the Bychkov-Rashba 



In what follows we consider the case when SO is small 
as compared to the magnetic splitting, so we can expand 
the self energy in a small parameter gkp/M <C 1. Then, 
we find 



Sf^k)-^ 2 1 + 



d 3 k' 



9 2 k 2 z \ 

AM 2 ) J (2tt) 



g' 2 k' 2 



AM 2 J £- £l (k')+M + ^ 



+ NiVr 



d 3 k' g' 2 k' 2 



1 



(2tt) 3 AM 2 e 



e 2 (k') + fj, + iS ' 
(56) 



and a similar expression for T, R {e,k). In these formulas 
we can put e = 0, k = kpitf), and average over the Fermi 
surfaces. This leads to the relaxation times for electrons 
corresponding to the bands 1 and 2, 



1 



n 



27TZ/1 Ni V 2 



1 



n 2 h 2 \ 

V i \2M 2 J ' 



(57) 



1 (, 

T2 = 2irv 2 N i V 2 y-^UM^J ' (58) 

where g-\ = g{kF*\) and g± = g{kn)- The above formulae 
explicitly show the contribution of SO interaction to the 
relaxation time. 

We also need to calculate the densities of states v\ and 
V2 at the Fermi surfaces and to relate them to the ones 
in the absence of the spin-orbit interaction, and v^. 
To do this we write 



/d 3 k 
j—y; 5 (n - e 1(2) {k)) 



and then, upon integrating over k x , we obtain 



^1(2) 



where 



_ (2m) 1 / 2 r 
(2^) 3 J e 



1(2) 



(k)<M (M-4( 2) (k)) 



1/2 



2m 



(59) 



(60) 



(61) 



Sk = 



{M-gk y y+g 2 k 2 



1/2 



(62) 



Expanding in powers of gkp/M -C 1 and calculating the 
integrals in (60), we find 



v 1 = v\ 1 



22 2 
m g^ m 9i 



9k 2 



2M 



(63) 



and 



v-i = v\ 1 + 



22 9 
m <7| mg| 



2k 



2M 



(64) 



The relaxation times T\ and t 2 enter the expressions 
for the renormalized retarded {G R ) and advanced {G A ) 
Green functions 



G RA {e,k) 



e — k 2 /2m + fi — Ma z + gk y a z — gk z a y 



(e - £ik + M ± */ 2r i) (e - e 2 k + M ± V 2t 2) 

(65) 

Now we have already all quantities necessary to calculate 
the renormalized vertex r aj 3ys. The general ladder- type 
equation for such a vertex reads 



+ N iVo n-auPs{u,q)r uslS {uj 7 q) 



(66) 
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where Vo is the matrix element of the short-range impu- 
rity potential, and 

_G« ( w ,k + q) Gl(0, -k) . 

(67) 



As before, we will restrict ourselves by considering the 
coupling constant g to be small as compared to the spin 
splitting, gkF/M < 1. 

Using Eqs. (66) and (67), and integrating over k x , we 
find for q, co = 



i(2m) 1/2 

iW»(o,o) = J 2 x 



dkydk- 



[s k - Ma 2 + .g(fc y cr z - k z a y )] av [s k -Ma z - g{k y a z - k z a y )] f 



dkydk- 

'4(k)< AI 



e' 1 ( k )<A' s k (s k + s_ k )(s-k - s k + i/ T i) [A* - + k 2 )/2m - s k 

[-s k - Mcr z + ff(fcj,cr z - k z a y )\ av [-s k - Mct z - .9(fc y (T2 



1/2 



^z CT y)]/3s 



Sk(sk + s_k)(sk - s-k + V r 2) [m - (^y + k 2 )/2m + s k ] 



1/2 



(68) 



If we put g = in Eq. (68), we find that only two matrix 
elements are nonzero, namely II -p -p -p -p (0, 0) and 11^1 (0,0). 
Assuming Dq 2 r, wr < 1, we can just supplement this 
result by the terms of the expansion in small q and lo 

II||f| (w, q) = 27r^Tj (l — D^q 2 T^ + iwr-f) , (69) 

and a similar one for <j). 

Let us consider now the limit of weak SO coupling 
and expand H au p s (q 7 uj) in powers of gkp/M. The first 
non- vanishing term of this expansion is quadratic in this 
parameter. According to Eq. (68), the other small pa- 
rameter in the limit of g is gkpr. Thus, if we assume 
Mr <C 1, wc can neglect the terms ~ (gkpr) 2 and keep 
the terms ~ gkp/M. In the opposite limit, Mr 3> 1 
(clean ferromagnet), we should keep terms ~ (gkpr) 2 
and neglect terms ~ gkp/M. Since our considerations of 
the localization corrections arc limited by epr 3> 1, the 
dirty case is possible only for some weak ferromagnets, 
when M <C ep (in other words, for ferromagnets with 
very low polarization). 



In view of Eq.(49), gkp is the amplitude of the spin- 
flip process. In the classical picture gkp is the angle of 
spin rotation in the unit time. Hence, gkp r is just the 
angle of the classical spin rotation due to the Bychkov- 
Rashba perturbation at the mean free path I of the elec- 
tron. Thus, the smallncss of the parameter gkp r corre- 
sponds to a small spin rotation angle at the length I. 

Since considerations in this paper are restricted to the 
case Mr 3> 1, wc will not discuss here the opposite limit, 
Mr <C 1. Such a case needs a special treatment. Indeed, 
as we pointed out in Sec. 1, only the condition Mr ^> 1 
allows to restrict oneselves to the triplet Cooper chan- 
nel. The result on the only nonvanishing matrix elements 
n TTTT (0,0) and n uu (0,0) for g = (see Eqs. (68) and 
(69)) also refers to the case Mr 3> 1 since in Eqs. (65) 
and (68) we account for the self energy as a small shift 
of poles of the Green function from the real axis. 

In the limit of small SO interaction, gkpr <C 1, and for 
Mr 3> 1, we can neglect in Eq. (68) all terms of the order 
of gkp/M. As a result, we find only two nonvanishing 
matrix elements 



IIf|||(w, q) = ITp (co», q) = 2-kv^t\ I 1 — Df q T| 



^ 2 2 ; 2 



(70) 



n llll( w ^) = n x (w,(7) = 27r^rj. ^1 - D i q 2 T l +iujT l - ^gfrf k 



2 Jl 7„2 



which give rise to the Cooperon r T (u;, q) in the form of Eq. (18), with the effective spin-orbit relaxation time 

1 4 



r so t 



■ 5| kpi Tj . 



(71) 



(72) 



Similar expressions can be obtained for r^(u;, q) and f so ^. 

If we take the limit gkpr » 1, then for Mr > 1 we still can consider the SO interaction as a small perturbation, 
gkp/M <C 1. Calculating the integral (68) up to the second order in this parameter, we find 



n ai//3s (w,g) 



37T^Ti 

8(.gifcF|Ti) 2 



1 — Diq 2 Ti + iujTi 



2M 3M 2 



m 2 g 2 
2k 2 Fl 



(1 - Cz)ai/(1 - 0- z )f3 s 
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9 _T£l a v a y 

3M 2 au 08 



+ 



9j k Fj 

6M 2 



5 av {l — <T z )f3s + 



9j k Fj 

6M 2 



(1 - o z ) av 8 i 



(3s 



mg\ 
M 



(1 - CT z )au(l + ^z)(3s 



8(g T fc F TT 2 ) 2 



1 - D 



t 2 + ioJT 2 + 



mg 2 



2M 3M 2 



m 2 g 2 
2k F ^ 



(1 + cr 2 )« y (l + Vz)ps 



3M 2 a a v°p s + 



6M 2 



)(3s 



9} ™F1 

6M 2 



(1 + Vz)av$0s + 



mg 2 



(1 



T z )f3s 



(73) 



Using (73), we can find all nonvanishing matrix elements of H av p s {uj, q) and then solve the ladder equations of Fig. 7. 
We do not make this calculation since we notice that the presence of the small factor l/(gk F T) 2 leads to complete 
suppresion of the Cooperon. Thus, the localization corrections exist only for gk F r < 1. 

Finally, we present the calculation of a classical correction to the conductivity due to the SO interaction in Hamil- 
tomian (49). After calculating the loop diagram with the Green functions (65), and using (57), (58), (63), (64), we 
find 



r(°) 



l T'T 



2 2 



m<? T V[ 9\ k Fl \ e 2 n in 



M 



v } 12M 2 



2 2 

m z g^ 

T 2 
K Fi 



m 9\ V] 9\ k F\ 



M 



V[ 12M 2 



(74) 



This expression contains corrections related to renormal- 
ization of the density of states (second and third terms in 
each bracket) and of the scattering time (fourth terms) . 

We should also take into account a renormalization of 
the chemical potential ijl. For this purpose we calculate 
the total number of spin up and down electrons and im- 
pose the condition of constant particle number. After 
simple calculations we find the correction to the chemi- 
cal potential 



A/1 = 



g^k F ^+glk 3 Fl m(g^k Fr +gfk F i) 



(75) 



6M(k F f + k Fi ) 2(k F] +k Fl ) 

Calculating the corresponding corrections to Fermi mo- 
menta, we find the corrections to spin up and down 
classical conductivities connected with renormalization of 
the chemical potential 

Act^ _ m(g 2 k 3 F1 + g 2 k 3 Fl ) 3m 2 (g 2 fc FT + g 2 k Fl ) 



V ZZ 



2Mk 2 F1 (k F1 + k Fl ) 2fc 2 T (fc FT + k Fl ) 



(76) 



Aal z \ _ m(g 2 k 3 n + g\ k Fl ) _ 3m 2 (g 2 . k F] + g 2 k F[ ) 
2Mk Fl (k F ^ + k F i) 2k F , (k F 

t + fc n) 



J® 

<Jzz 



(77) 

It should be noted that the classical corrections, (76) and 
(77), can be of any sign since the inequality M 3> 1/r 
does not necessarily imply any relations between the 
magnitude of M and the Fermi energies of majority and 
minority electrons. 

The total correction to conductivity includes all cor- 
rections to the classical part, Eqs. (74), (76) and (77), as 
well as the quantum correction in the form of Eqs. (21) 
or (36 a,b) with f so , given by (72) for gk F r < 1. 



VII. EFFECT OF INTERNAL MAGNETIC 
INDUCTION IN A FERROMAGNET 



In our model we have taken into account the effect of 
magnetic induction B, which enters the kinetic energy 
of electrons through a vector potential A and leads to 
the suppression of localization corrections, Eqs. (22) and 
(36b). The total magnetic induction B inside a ferro- 
magnet includes the external magnetic field H ext and 
the internal magnetic induction B m * , B = H ext + B m * , 
where B™ 1 = Air (Sij — n,j ) Mqj , Mo is the magnetiza^. 
tion vector and ntj is the demagnetizing factor tensor.EEl 
The magnitude of the internal magnetic induction may be 
rather high in strong ferromagnets. Nevertheless, weak 
localization corrections were observed, e.g., in Ni films.E^I 

We can present some numerical estimations of the mag=. 
nitude of B mt . For example, in the case of Fe we takeEa 
4ttM ~ 2 T. Thus, for the bulk Fe, when the demag- 
netizing field is negligibly small, we have B mt ~ 2 T. 
The critical magnetic induction B crlt , which can totally 
suppress the localizations correction, is determined dif- 
ferently in 3D, and 2D cases. In the 3D case we can esti- 
mated it byM Ib — I (I is the electron mean free path). 
For the parametecs^of very pure bulk Fe [m ~ 4 mo and 
t ~ 5 x 10~ 13 s),E3 we find I ~ 4 x 10~ 5 cm, and con- 
sequently B cr n ~ 50 Oe. This estimation shows that 
the localization corrections in pure bulk Fe are totally 
suppressed by the internal magnetic induction. On the 
other hand, in not so pure metals, or in magnetic al- 
loys, and amorphous materials, one can expect a much 
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shorter mean free path, like I ~ 10 -6 cm, which is still 
large enough for the localization corrections to be small 
(kpl ;§> 1). In case of I = 1CP 6 cm, one finds B cr n — 7 T, 
which is significantly larger than the internal induction 
of 2 T and, therefore, makes the localization corrections 
observable. 

In the case of thin magnetic films, the demagnetizing 
factor is of crucial importance. For example, when the 
magnetization vector Mo is perpendicular to the inter- 
face, the demagnetizing factor is unity, n zz — 1 (z axis is 
perpedicular to the plane; other components of demag- 
netizing tensor are very small), and we have B m * = 0. 
In this case, the results of Sec. 4 can be applied with 
the magnetic induction B equal to the external magnetic 
field H ext . The critical magnetic field, oriented perpen- 
dicular to the plane, which suppresses the localization 
correction, can be estimated as in the 3D case. 



On the other hand, in the case of in-plane magneti- 
zation, the demagnetizing factor is much smaller than 
unity, and the non- vanishing internal magnetic induction 
B mt = 47rMo is directed along the film. However, the ef- 
fect of parallel magnetic induction on the localization cor- 
rections in a strongly two-dimensional system is absent 
(Sec. 4) . For a quasi-2D system like a quantum well, the 
effect of parallel induction is non- vanishing but weak. We 
can write the corresponding expressions for localization 
corrections in the presence of in-plane magnetic induction 
by simply generalizing the results of Refs. [ |36| , [37| . 

If the film thickness L is large with respect to 
the electron mean free path but still -small enough 
to consider the film as a 2D system|m I L 

Hpfso)^ 1 + (Dt^)^ 1 ] , we can obtain, using Ref. [ 
Bqj, the dependence of localization correction on mag- 
netic induction in the limit of small B, i.e., for l B 3> L: 



Aa(B) MO) = ^ 



L 2 Di 



In 1 



so] 



L 2 D 



121% 



In 1 



(>1 



+ T 



121% 



(78) 



It shows a weaker dependence on the magnetic induction, as compare to the 3D case. For larger values of B, when 
Ib < L, the dependence on magnetic induction is as for the 3D case, Eq. (22). The critical value B cr i t can also be 
estimated as for the 3D case. Thus, for thick clean magnetic films, such that L ^> I, and l> Iq= [c/ (AckMo)] 1 / 2 , the 
in-plane internal magnetic induction suppresses the localization corrections completely. But in dirty or amorphous 
thick films with I < lo, they can be observed. Using the parameters of Fe, we find Iq — 1-8 X 10~ 6 cm. 

If the film thickness L is smaller than the mean free path, I (ballistic regime) , we use the result of Ref. [ |3^] , which 
can be presented in a simple form for some intervals of BEB 
To avoid combersome formulae, we introduce the following notations 



1 



1 



1 



T cT,l r s°T,i T vt,l 



L c1,l = (• D T4 r cT,i) 1 2 • 



(79) 



When L c , > !, we obtain 



Aa(B) - Aa(0) = — 2 



161% r T 



161% n 



(80) 



When L c , I > l 2 B /L, 



Aa(B) - Acr(0) 



4tt 2 



3/| n 



(81) 



and when I, l B /L ^> L c , 

' . <82) 

We also find that the critical in-plane magnetic induc- 
tion in the case of L <C I can be estimated from l B ~ L, 
which gives rise to a much larger critical value of mag- 
netic induction, Bg^jB?" ~ {l/Lf » 1. In other 



words, even for clean magnetic films, the in-plane mag- 
netic induction does not suppress completely the local- 
ization corrections, if the film width is sufficiently small. 

In view of the above considerations, the best configu- 
ration for observing weak localization effects in a ferro- 
magnet is to use a thin film with perpendicular easy axis, 
and apply a perpendicular magnetic field. 
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VIII. SUMMARY AND CONCLUSIONS 

We have analyzed the localization corrections to elec- 
trical conductivity in magnetically polarized materials 
with spin-orbit interactions. The strong magnetic po- 
larization excludes processes with the singlet Cooperon, 
which are responsible for the antilocalization in nonmag- 
netic materials with SO scattering. As a result, the quan- 
tum correction to conductivity is always negative in fer- 
romagnets and leads to negative magnetoresistance. 

The strength of SO interaction, together with the 
phase relaxation time due to inelastic processes, de- 
termine the magnitude of these corrections. In the 
case when the SO interaction is associated with scat- 
tering from impurities and/or other defects, the effec- 
tive SO scattering time f so entering the Cooperon, de- 
pends on the dimensionality of the system and on the 
magnetization orientation with respect to the plane of 
the system (in the case of two-dimensional or quasi- 
two-dimensional systems). In the case of strongly two- 
dimensional ferromagnets with in-plane magnetization, 
the inverse time l/f so is zero. This increases the magni- 
tude of the localization correction. The vanishing value 
of of 1 /f so is essentially related to the spin-flip scattering 
of the Cooperon; the usual contribution from the spin- 
conserving scattering is cancelled by the spin-flip contri- 
bution, which enters the Cooperon with an opposite sign. 
In the quasi-two-dimensional case both contributions are 
present, but they do not cancel each other. 

We have also found the effective spin relaxation time in 
the case of Bychkov-Rashba SO interaction. It contains 
contributions from both spin-flip and spin-conserving 
scattering processes. 

We think that good candidates for observations of the 
localization corrections— are also semiconducting ferro- 
magnets like GaMnAsjla which are recently extensively 
studied in view of their possible applications in spintron- 
ics. Another example is a new ferromagnetic semicon- 
ductor CaBg, where a very small magnetization can be 
combined with a small electron density.cJ 
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